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ABSTRACT 

We study the distribution of interstellar matter near the galactic plane on the basis of open 
star clusters and reviewed the correlation of extinction with different physical parameters 
of the clusters to understand the galactic structure. It is seen from the extinction towards 
open clusters that about 90% of the absorbing material lie within —5° < b < 5° of the 
galactic plane. The mean thickness of the absorbing material, which is determined in terms 
of half-width value /3, is estimated to be about 125 ± 21 pc. We show that the interstellar 
absorption follows a sinusoidal variation with galactic longitude and maximum and minimum 
absorptions occur at I ~ 48° ± 4° and I ~ 228° ± 4° respectively. It is found that the galactic 
plane defined by the reddening material is inclined by an angle of 0°.6 ± 0°.4 to the formal 
galactic plane and inclination is maximum at I ~ 54° ± 6°. The reddening analysis has been 
used to constrain the Solar offset which is found to be about 22.8 ±3.3 pc above the reddening 
plane. We obtained a scale height of 53 ± 5 pc for the distribution of open clusters while it is 
186 ± 25 pc for the distribution of reddening material from the reddening plane. 

Key words: ISM: dust, extinction - Galaxy: structure, open clusters - method: statistical 



1 INTRODUCTION 



Since long time, distribution of interstellar material has been used 
to understand the galactic structure. A number of studies have been 
carried out on the basis of reddening variation of stars and various 
parameters like reddening plane, Solar offset, distribution of inter- 
stellar material in the sky etc. were derived using different kind of 
objects. Fernie (1968) studied the galactic structure on the basis 
of classical Cepheids. FitzGerald (1968) analysed the interstellar 
extinction within a few kpc of the Sun using colour excesses of 
7835 O to M stars and found that the reddening material is highly 
concentrated to the galactic plane. Neckel & Klare (1980) derived 
extinctions and distances for more than 1 1000 O to F stars and in- 
vestigated the spatial distribution of interstellar extinction in the 
direction of |6| < 7° .6. Lately, open star clusters have been used to 
study the distribution of interstellar extinction at low-galactic lat- 
itude rather than the stars. This kind of study has an advantage 
over the stellar study since the distances and reddening derived 
through open clusters are regarded as more reliable than those de- 
rived through individual stars. Furthermore the open star clusters 
span a wide range of galactic longitudes, latitudes and distances 
from the Sun. Janes & Adler (1982) studied the galactic structure 
and distribution of clusters in the galactic plane from a catalog of 
434 open clusters. Pandey & Mahra (1987, hereafter PM87) used 
photometric data of the 462 open clusters within 1.5 kpc of the Sun 
to investigate the spatial distribution of interstellar extinction and 
concluded that the galactic plane defined by the reddening material 
is inclined by an angle of 0°.8 ± 0°.2 and the Sun is situated at a 
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distance of about 10 pc above the plane of symmetry defined by the 
interstellar matter. Pandey et al. (1988), on the basis of young open 
clusters of the same sample, found these values as 0°.5 ± 0°.4 and 
28 ± 5 pc respectively. From the same data, though, Janes et al. 
(1988) concluded that there are no significant undulations around 
the plane of symmetry within 4 kpc of the Sun. 

However, conclusions drawn from the photometric data of the 
open clusters may have been affected by the considerable uncer- 
tainty in their parameter determinations like distance, reddening, 
age etc. as well as lack of homogeneous data. During last ten years, 
new data on the open clusters along with more precise determi- 
nation of physical parameters have become available mainly due 
to large size telescopes and improved detectors like CCD which 
allowed us to obtain deeper photometry of the clusters. Chen et 
al. (1998) derived an analytical expression for the interstellar ex- 
tinction as a function of galactic longitude and distance based on 
the observed distances and colour excesses of 434 open clusters 
within \b\ < 10°. Recently Dias et al. (2002) provided a com- 
prehensive catalogue of more than 1600 open clusters which also 
include the information available in the Lynga catalogue (1987) 
as well as WEBDA catalogue (Mermilloid, 1995). The catalogue 
contains almost all the clusters known till date. It therefore seems 
worthwhile to re-investigate the classical problem of galactic struc- 
ture with this improved informations. In the present paper, we in- 
vestigate the spatial distribution of interstellar material and study 
the Galactic structure using the physical parameters derived for a 
large sample of open clusters. The detail about the data used is ex- 
plained in Sect. 2. In Sect. 3, we analyse the spatial distribution of 
interstellar extinction while the Sect. 4 deals with the correlation of 
interstellar extinction with the different physical parameters of the 
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Figure 1. The distribution of open clusters in the I — b plane is shown in 
(a). A number distribution of the clusters in log(Age), distance from the Sun 
(cumulative number) and the distance perpendicular to the galactic plane are 
shown in (b), (c) and (d) respectively. Few odd clusters lying quite far-away 
in the x-axis are not shown in (d) for the clarity of the histogram. 

clusters in order to understand galactic structure. The main results 
of the present study are summarized in Sect. 5. 
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Figure 2. The distribution of open clusters as projected on the galactic 
plane. Points are normalized to the extinction values such that the smallest 
and largest points denote E(B — V) = 0.01 mag and 2.55 mag respec- 
tively. The position of Sun is shown by a symbol of star at the center. Five 
concentric circles at an equal distance of 1 kpc are drawn. The number of 
clusters in each distance bin are also given. 

rial may be lying below the galactic plane. A further study of this 
effect shall be discussed in Sect. 4.2. 



2 THE DATA 

Data for the present work is taken from the data base of open star 
clusters compiled by W. S. Dias (Dias et al. 2002) at the web site 
http://www.astro.iag.usp.br/~wilton/ Dias et al. (2002) reported a 
catalogue of 1599 clusters at the CDS, however, a total of 1632 
clusters are available in the online catalogue updated till March, 
2005 in which 747 clusters have distance and reddening informa- 
tions. Out of them, we have included only 722 clusters for our anal- 
ysis since 20 of them are considered to be doubtful by the DSS 
images inspection (Dias et al. 2002) while 5 are POCR (Possible 
Open Cluster Remnant, Bica et al. 2001). The age informations of 
602 clusters are also available out of these 722 clusters. The clusters 
are almost uniformly distributed in the longitude / but are densely 
packed in the latitude b near the galactic plane as seen in Fig. 1(a). 
The age histogram drawn in Fig. 1(b) shows that the clusters are 
known in all the age range from 6 Myrs to as much as 10 Gyrs. 
Though the clusters are observed up to a distance of about 15 kpc, 
however, it is seen that the number of clusters decreases with dis- 
tance which suggests of possible distance and/or reddening effect 
on the detectability of clusters. When we plot the cumulative dis- 
tribution of the clusters as a function of distance from the Sun (see 
Fig. 1(c)), we found that the present cluster sample is complete only 
up to a distance of about 2 to 3 kpc. It is seen from the histogram of 
the distance perpendicular to the galactic plane (z = r sin&) plot- 
ted in Fig 1(d) that almost 90% clusters lie within 200 pc of the 
galactic plane. We notice that the clusters distribution is maximum 
at z ~ —25 pc which indicates that the maximum reddening mate- 



3 INTERSTELLAR EXTINCTION TOWARDS OPEN 
CLUSTERS 

One of the characteristics of the interstellar extinction is its irreg- 
ular structure in the sky. Using the reddening of 722 open clusters 
which have been determined independently through the analysis of 
colour and magnitude of stars in the clusters, we map the extinc- 
tion variation in the sky that is shown in Fig. 2. All the clusters are 
put on the galactic plane to draw the map. Since most of the clus- 
ters having reddening information are observed within 5 kpc of the 
Sun, we map here only first 5 kpc regions shown by 5 concentric 
circles at an equal distance of 1 kpc around the Sun. The size of 
the dots are an indication of extinction in the sense that bigger the 
dots, larger the extinction is. The map shows a complex structure of 
the sky. One of the striking feature of the map is clumpy behaviour 
of the cluster distribution at some places like I ~ 320° — 350°. 
The concentration of the reddening material to the galactic plane 
is quite visible in the region I ~ 330° .. 0° .. 150° while the least 
absorption lie in the region I ~ 180° to 270°. A galactic dust dis- 
tribution map up to 3 kpc distance from the Sun shown by Neckel 
& Klare (1980) also show similar picture. 

We compared the E(B — V) values for the open clusters given 
in the Dias et al. (2002) with the extinction map given by Schlegel 
et al. (1998, hereafter SFD98) which is based on the IRAS 100 
fim surface brightness converted to extinction and a comparison 
in three different latitude bins are shown in Fig. 3. We found a 
large discrepancy in the extinctions values towards open clusters 
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Figure 3. A comparison of the interstellar extinction towards opon clusters 
with the Schlegel map for 3 different latitude range. The dashed lines denote 
the y = x line. 
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Figure 4. The variation of E(B — V) vs I (left) and b (right). In the upper 
panels, the histogram of the mean value of reddening using bin sizes of 10° 
in I and 1° in b are shown. 



in the SFD98 reddening map below b < 5°. It is not unexpected 
since most contaminating sources have not been removed from 
the SFD98 maps as well as temperature structure of the Galaxy 
is not well resolved at low galactic latitude. The comparisons at 
b > 5° show that the SFD98 map is overestimated by a factor of 
1.04 ± 0.09 and 1.12 ± 0.05 in comparison to the extinction based 
on the open clusters in the regions 5° < b < 10° and b > 10° re- 
spectively. Chen et al. (1999) also found that SFD98 reddening map 
overestimates E(B — V) by a factor of 1.16 and 1.18 based on the 
study of globular and open clusters respectively at |&| > 2.5°. In 
recent times, various authors have compared interstellar extinction 
in different regions of the sky with the SFD98 reddening map (cf. 
Arce & Goodman, 1999, Cambresy et al. 2001, 2005, Dutra et al. 
2002, 2003a, 2003b) and supported the fact that SFD98 map over- 
estimates the extinction in many parts of the sky, if not in the whole 
sky. 

In order to study the distribution of extinction along the galac- 
tic longitude and latitude, we plot E(B-V) against / and b in 
Figs. 4(a) and 4(b); and distribution of the mean value of E(B — V) 
in the 10° bin in I and 1° bin in 6 in Figs. 4(c) and 4(d) respectively, 
however, around the peak in Fig. 4(d), we made smaller bin size of 
0.5° due to higher concentration of clusters. We have calculated the 
mean value of reddening (E(B — V)), concentration of reddening 
(Cred) which is defined as E(B — V) in the sub-region divided by 
total E(B — V) of the clusters and number of clusters (Ncl) for 
5 different latitude range which are given in Table 1. From the fig- 
ures and table, one can easily notice following approximate though 
important features: 

(i) That the reddening is comparatively larger and scattered in 
the direction of galactic center than the small and systematic varia- 



Table 1. The distribution of the clusters and reddening in different latitude 
range. 



range in |6| 


E(B - V) 


Cred 


N CL 


(deg) 


(mag) 


(%) 




0-2 


0.59 


67 


383 


2-5 


0.40 


22 


190 


5-10 


0.29 


7 


84 


10-20 


0.19 


3 


54 


>20 


0.08 


1 


11 



tion in anti-center direction (see Figs 4(a) and 4(c)). On an average, 
the reddening is almost twice in the region I ~ 330° .. 0° .. 150° 
than the rest of the region. 

(ii) Most of the open clusters are concentrated near the galactic 
plane. However, we notice that the peak of the reddening is shifted 
by ~ 0.5° in the southwards of the galactic plane (see Figs. 4(b) 
and 4(d)). 

(iii) A simple approximation on the basis of E(B — V) infor- 
mations of the open clusters indicate that about 90% reddening ma- 
terial lie within b ± 5° of the galactic plane (see Table 1). 

It is seen from the age distribution of clusters under study that while 
young open clusters are confined very close to the galactic plane 
and have a clumpy distribution, the old age clusters are mostly 
distributed in the far-off region from the galactic plane. Our study 
shows that the reddening layer is thicker near the galactic plane and 
about 80% of clusters sample and about 90% reddening material lie 
within —5° < b < 5° of the galactic plane. We, therefore, consider 
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Figure 5. The normalized extinction E„(B — V) versus distance r for 36 zones of galactic longitude. Range of the longitudes for each zone are shown at the 
top of each plot. 



only those clusters for the further analysis which are lying in this 
region of the Galaxy. 

4 ANALYSIS AND RESULTS 

4.1 Variation of extinction with distance 

It is a well known fact that the interstellar matter is extremely in- 
homogeneous in the sky, therefore, we study the behaviour of ex- 
tinction variation in the sky by dividing it in small zones. In some 
previous studies of the galactic structure which have been carried 
out on the basis of extinction variation in the stars (e.g. Neckel & 
Klare, 1980), the sky has been split up in very small sub-regions. 
Chen et al. (1998), on the way to construct extinction map in the 
galactic plane below |fe| < 10° divided the cluster sample into 36 
cells with Al — 10°. Owing to the limitations of smaller sample 
of the clusters than the stars, we also divide our sample of clusters 
located in \b\ < 5° into 36 zones of 10° interval in longitude. We 
normalize the extinction by distance and a term, E„ (B — V) , called 
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Figure 6. The mean value of the normalized extinction as a function of I. 
A best fit sinusoidal variation is shown by continuous line. 
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Figure 7. The variation of as a function of z for 8 different regions in longitude which are written at the top of each diagram. A least square Gaussian fit 
around the maximum absorption is also drawn. 



normalized extinction has been evaluated as 

r 

where r is the distance of cluster from the Sun. A variation of 
E n (B — V) with the distance from the Sun for all the 36 zones 
are drawn in Fig. 5. In general, we see following features in the 
diagrams: 

(i) Within individual zones, the reddening shows a considerable 
but irregular variation which cannot be approximated by any simple 
analytical function. 

(ii) The scattering in extinction variation is more in the direc- 
tion of galactic center where most of the complex dust clouds are 
believed to be located. 

(iii) The E n (B — V) changes significantly from one zone to 
another zone from about 0.1 to 0.6. A similar variation has also 
been noticed by Chen et al. (1998) near the galactic plane. 

(iv) There is a general trend in the variation of E n (B — V) that 
it is maximum around the longitude range 30° — 50° and minimum 
around the longitude range 220° — 250°. This kind of variations 
are also seen by Arenou et al. (1992) and Chen et al. (1998) in their 
extinction models. 

We determine the mean value of normalized extinction for each 
zone which is plotted as a function of galactic longitude in Fig. 6. 



The variation can be described as a sinusoidal function and a least 
square fit gives 

E n (B- V) = 0.33 + 0.16 sin(Z + 43) 
±0.01 ±0.02 ±6 

where E n (B — V) is in mag/kpc and I in degree. From the equation, 
it is estimated that the extinction varies rapidly with distance in 
the direction of I ~ 50° while it varies with a very slow rate in 
the direction of I ~ 230°. The mean extinction in the region of 
160° < I < 300° is only ~ 0.07 mag/kpc which suggests that the 
sky is relatively transparent in this direction of the Galaxy. 



4.2 Variation of absorption with z 

To study the distribution of interstellar extinction perpendicular to 
the galactic plane, we initially determined the interstellar absorp- 
tion k as 



r 



R 



E(B - V) 



where R is the total-to-selective absorption which we have taken as 
3.1. Therefore, above equation can be expressed as 

k = Z.lE n {B - V) 
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Figure 8. Maximum absorption fco as a function of I. A least square sinu- 
soidal fit is drawn by a continuous line. 



Since the different parts of the sky do not behave similarly, we 
therefore divide the sky in 8 zones of the galactic longitude. The 
boundaries of these zones are chosen in such a way that the larger 
boundaries are made towards galactic center where smaller num- 
ber of clusters are seen per longitude bin due to higher extinction 
while smaller boundaries are made towards galactic anti-center due 
to opposite behaviour. Though it makes irregular boundaries but 
render lesser scatter within the fields. The boundaries of the re- 
spective zones are given at the top of each panel of Fig. 7 where 
we show the variation of absorption in all the 8 zones as a function 
of z. Here we have plotted the running mean of fc for 60 pc inter- 
val of z against the mean value of z of that interval. A Gaussian 
fit around the central peak value of z in all the diagrams of Fig. 7 
is also drawn. Although the distribution of interstellar absorption 
varies significantly from one zone to another but the strong con- 
centration of absorption in the galactic plane is clearly visible in all 
the zones except 230° < I < 300°. Guarinos (1991) has pointed 
out that between Orion and Sco-Cen complex of the above zone, 
the interstellar medium is transparent around the Sun up to at least 
1.0 to 1.4 kpc. Neckel & Klare (1980) reported lowest extinction in 
the region 210° < I < 250° from the analysis of early type stars 
within 1 kpc of the Sun. Chen et al. (1999) also noticed deviation 
in the region 225° < / < 290° in order to constrain the extinction 
model based on the COBE/IRAS all sky reddening map. 

The thickness of the absorbing layer in different zones can be 
expressed in terms of half-width value, (3, which is defined as the 
separation of z values at the 1 /e of the maximum value of inter- 
stellar absorption (Neckel, 1966). We determined thickness for the 
each zone and a final value of the (3 has been derived as a mean 
of the thickness which is estimated to be 125 ± 21 pc in the z di- 
rection. FitzGerald (1968) has reported /3 in the range of 40 to 100 
pc on the basis of colour excess versus distance diagram of about 
8000 O to M type stars while studying open clusters in |&| < 10°, 
PM87 found it greater than 100 pc for most of the sub-regions of 
different longitude range in the sky. 

In order to analyse the absorption variation in more detail, we 
have determined the maximum value of absorption fco and corre- 
sponding distance from the galactic plane zq from the fit. The vari- 
ation of fco as a function of mean value of longitude for each zone 
is drawn in Fig. 8 which shows a sinusoidal kind of variation. A 
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longitude (deg) 

Figure 9. The height z above or below the galactic plane at the maximum 
absorption fco as a function of /. A least square sinusoidal fit is shown by a 
continuous line while dashed line shows the mean value. The lone point in 
the vicinity of I ~ 250° drawn by open circle is omitted from the fit. 

least square solution to the sinusoidal function gives 

fco = 1.30 + 0.69 sin(/ + 42) 
±0.03 ±0.04 ±4 

The absorption is found to be maximum in the direction of / ~ 
48° ± 4° and it is minimum towards I ~ 228° ± 4°. Fernie (1968) 
has also found a maximum average absorption towards I ~ 50° 
based on the study of classical Cepheids. Based on the study of 
open clusters, Lynga (1982) and PM87 have found I ~ 230° and 
/ ~ 50° for the minimum and maximum absorptions respectively. 
Using an inverse technique of Tarantola & Valette (1982) on open 
clusters, Chen et al. (1998) found the lowest absorption in the di- 
rection of I = 210° to 240°, however, they reported the highest 
absorption towards I = 20° to 40° which is slightly lower than 
those found in all other studies. 

To study the variation of maximum absorption perpendicular 
to the galactic plane, we plotted the value of zo against the mean 
value of longitude in Fig. 9. A least square solution for a sinusoidal 
function gives 

z = -22.8 + 41.1 sin(2 + 36) 
±3.3 ±4.8 ±6 

We have excluded the lone point around 230° < I < 265° in 
drawing the fit where zo has quite a large shift from the sinusoidal 
function. From the figure, we draw following conclusions: 

(i) That the maximum extinction lies above the galactic plane to- 
wards galactic center and below the galactic plane towards galactic 
anti-center. 

(ii) It shows that the distance of the galactic plane at maximum 
absorption is symmetric around z ~ —22.8 pc which indicates that 
the larger amount of reddening material is lying below the galactic 
plane which is in confirmation with our approximate estimation in 
the earlier section. 

(iii) A symmetric variation around z ~ —22.8 pc means that 
the galactic plane of symmetry defined by the reddening material is 
lying below by the same value to the formal galactic plane (b = 0° 
plane). This suggests that the Sun probably lies about 22.8 ± 3.3 
pc above the galactic plane defined by the reddening material. The 
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Table 2. The value of Solar offset above the reddening plane in previous 
studies. 



Zq (pc) 


Based on the study 


Reference 


30 


15 ±3 


Wolf-Rayet stars 


Conti & Vacca (1990) 




15.5 ±0.7 


IRAS point-source counts 


Martin (1995) 







20.5 ±3.5 


optical star counts 


Humphreys & Larsen (1995) 


ft 


~ 14 


near-IR surface brightness map 


Binney et al. (1997) 


N 

-30 


10 - 12 


OB stars within \b\ < 10° 


Reed (1997) 




27 ±3 


star catalogue and simulation 


Mendez & van Altena (1998) 




27.5 ±6 


extinction model 


Chenetal. (1999) 


60 


24.2 ± 1.7 


0-B5 stars 


Maiz-Apellaniz (2001) 




22.8 ±3.3 


open clusters within 6 < 5° 


present study 


-90 



offset of the reddening plane from the galactic plane has impor- 
tant bearing on the determination of density distribution of different 
kind of stars, particularly young stars, which are closely situated 
near the galactic plane (Mendez & van Altena 1998). Numerous 
studies have been carried out to determine Solar offset from the 
reddening plane using different kind of objects (see Table 2 for few 
of them). Most of the recent studies show a shift in the range of 
15 — 30 pc in the north direction of the galactic plane. Our de- 
termination of Solar offset is thus in close agreement with these 
estimates. 

(iv) The galactic plane seems to be tilted and maximum upward 
tilt is found to be in the direction of / ~ 54° ±6° . PM87 and Pandey 
et al. (1988) have also reported an upward tilt in the direction of 
/ ~ 60° and 50° respectively. 



4.3 Reddening plane versus galactic plane 

It is a well known fact that the formal galactic plane does not co- 
incide with the planes of symmetry defined by the different objects 
like neutral hydrogen layer (Gum et al., 1960), Cepheid variables 
(Fernie, 1968), WR stars (Stenholm, 1975), HII regions and su- 
pernova remnants (Lockman, 1977), the molecular clouds (Cohen 
& Thaddeus, 1977, Magnani et al., 1985), Open clusters (Lynga, 
1982, Pandey et al. 1988) etc. To examine the relationship between 
reddening plane and formal galactic plane, we also investigate the 
distribution of z in a 1 kpc bin as a function of r following the ap- 
proach given by Fernie (1968). This we determined in two different 
directions having maximum absorption i.e. towards I = 54° ± 35° 
and I = 234° ± 35°. A variation in the mean value of z as a func- 
tion of distance is plotted in Fig. 10 where distances towards I = 
54° ± 35° are taken with positive sign while distances towards I = 
234° ± 35° are taken with negative sign. If zq is the distance of the 
Sun from the reddening plane and <j> is the angle between redden- 
ing plane and formal galactic plane then the relation between two 
planes can be expressed as 

2 = r sin (f> + zq 

where both z and r are expressed in pc. Since we have earlier de- 
termined 2 Q = 22.8 ± 3.3 pc, we put this value here and a least 
square regression fit gives us 

<j> = 0.6° ± 0.4° 

This shows that the reddening plane is inclined by an angle of 
0.6° ±0.4° with respect to the galactic plane. PM87 found an angle 
of 0°.8 ± 0° .2 between these two planes while Pandey et al. (1988) 
reported it as 0° .5 ± 0° .4. All the results are thus consistent within 
their quoted errors. If we assume that the Sun lies at a distance of 




Figure 10. Mean value of the 2 in a 1 kpc bin as a function of distance 
in the direction of I = 54° ± 35° (taken as a positive distance) and I = 
234° ± 35° (taken as a negative distance). The least square fit is shown by 
a continuous line. 



8.5 kpc from the galactic center then we estimate that the reddening 
plane may cut the formal galactic plane at a distance of 6.3 ± 0.9 
kpc from the galactic center. 



4.4 Determination of scale heights 

To determine the scale height of the distribution of open clusters 
from the reddening plane, we defined the plane of symmetry due to 
reddening material as 

z = 2 Q + r sin & cos — r cos b sin <f>cos(l — 9 t ) 

where <j> is the inclination angle between two planes and 8 t is the 
angle of maximum inclination of the reddening plane. Other sym- 
bols have their previous meanings. Suppose a large scale volume 
density of the reddening material, n(r), follows a decaying expo- 
nential law with respect to the reddening plane then it can be ex- 
pressed as 

n(z') = no exp( — \z'\/h) 

where h is the scale height of the open cluster distribution from the 
reddening plane. The mean value of z has been calculated in 20 pc 
intervals of z' and number of open clusters are counted in each bin. 
In Fig. 1 1, we have plotted the logarithm of number distribution of 
star clusters as a function of z . From a least square linear fit, we 
determined 

h = 53 ± 5 pc 

The percentage of clusters within z of the reddening plane can be 
determined as 

n(z') = [1 - exp{-\z'\/h)\ x 100% 

It shows that about 90% of the open clusters are confined within a 
distance of 122 ± 12 pc of the reddening plane. 

Similarly, to determine the scale height of the distribution of 
interstellar material normal to the reddening plane, we draw loga- 
rithm of the mean value of k in 20 pc interval of z' as a function of 
z' in Fig. 12. We further fit a linear regression line of the form 

k(z') = fco exp(— \z'\/h') 
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Figure 11. The density distribution of clusters as a function of z' . The 
vertical axis is in the natural logarithm of z'. Continuous line is a least 
square fit to the points. 
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Figure 12. The density distribution of absorption as a function of z' . The 
vertical axis is in the natural logarithm of z' . Continuous line is a least 
square fit to the points. 

where b! is the scale height of the interstellar material from the 
reddening plane. A least square fit gives 

ti = 186 ± 25 pc 

From the above equation, it is found that about 90% of the redden- 
ing material lie within a distance of 430 ± 60 pc from the reddening 
plane. 

Pandey et al. (1988) estimated a scale height for the distribu- 
tion of open clusters h — 52 ± 6 pc from their study of nearby 
open clusters while Lynga (1985) reported a scale height of 60 pc 
based on the study of young open clusters. Our value is, thus, in 
agreement with the earlier estimates. The scale height of the dis- 
tribution of reddening material determined by PM87 as 160 ± 20 
pc which is smaller than our result but their value is estimated with 
respect to the formal galactic plane which is shifted by a distance 
of 22.8 ± 3.3 upwards from the reddening plane. 



Table 3. Summary of the parameters derived in the present study 



Parameters 


Value 


Distance of the Sun above the reddening plane 


22.8 ± 3.3 pc 


Direction of maximum absorption 


I ~ 48° ± 4° 


Angle of inclination between reddening plane 


0°.6±0°.4 


and formal galactic plane 




Direction of maximum inclination 


I ~ 54° ± 6° 


Half- width value of the reddening material 


125±21 pc 


Scale height of the distribution of open clusters 


53±5 pc 


Scale height of the distribution of reddening material 


186±25pc 



5 SUMMARY 

In order to study the distribution of the interstellar extinction in 
the sky and investigate the galactic structure in this context, we 
used the data base of 722 open star clusters available at the web- 
site http://www.astro.iag.usp.br/~wilton/ which is complied by the 
Dias et al. (2002). On the basis of reddening of open clusters within 
|b| < 5°, we found that the extinction varies with distance in a very 
slow rate of about 0.07 mag/kpc in the region of 160° < / < 280°, 
however, it is highly variable in other parts of the Galaxy. We de- 
termined that almost 90% of the interstellar material lie within 5° 
of the galactic plane and the maximum extinction lies above the 
galactic plane towards galactic center and below the galactic plane 
towards anti-center direction. The mean thickness of the reddening 
material in the galactic plane, which is described in terms of the 
half-width value f3, is found to be 125 ± 21 pc. An analytical ex- 
pression for the interstellar extinction as a function of galactic lon- 
gitude and distance has been derived. The reddening analysis has 
been used to constrain the Solar offset from the galactic plane. We 
found that the galactic plane defined by the reddening material is 
symmetric around the z = —22.8 ± 3.3. It is seen that the galactic 
plane defined by the distribution of reddening material is inclined 
by an angle of 0°.6 ± 0°.4 to the formal galactic plane and the 
inclination is maximum at / ~ 54°. We summarize main results 
of our study in Table 3. Though most of our results are in close 
agreement with the values determined in previous studies but it can 
be stated that a larger sample of clusters with better knowledge of 
their parameters enable us to achieve a considerable improvement 
in the analysis of galactic structure. Furthermore, we provided an 
overview of the distribution of reddening material in the galactic 
plane. 
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